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Abstract 

We regularize compact and non-compact Abelian Chern-Simons-Max- 
well theories on a spatial lattice using the Hamiltonian formulation. We 
consider a doubled theory with gauge fields living on a lattice and its 
dual lattice. The Hilbert space of the theory is a product of local Hilbert 
spaces, each associated with a link and the corresponding dual link. The 
two electric held operators associated with the link-pair do not commute. 
In the non-compact case with gauge group M, each local Hilbert space 
is analogous to the one of a charged “particle” moving in the link-pair 
group space M 2 in a constant “magnetic” background held. In the com¬ 
pact case, the link-pair group space is a torus U( l) 2 threaded by k units 
of quantized “magnetic” hux, with k being the level of the Chern-Simons 
theory. The holonomies of the torus I/(l) 2 give rise to two self-adjoint 
extension parameters, which form two non-dynamical background lattice 
gauge helds that explicitly break the manifest gauge symmetry from U( 1) 
to Z (k). The local Hilbert space of a link-pair then decomposes into rep¬ 
resentations of a magnetic translation group. In the pure Chern-Simons 
limit of a large “photon” mass, this results in a Z( k )-symmetric variant 
of Kitaev’s toric code, self-adjointly extended by the two non-dynamical 
background lattice gauge helds. Electric charges on the original lattice 
and on the dual lattice obey mutually anyonic statistics with the statis¬ 
tics angle ^. Non-Abelian U{k) Berry gauge helds that arise from the 
self-adjoint extension parameters may be interesting in the context of 
quantum information processing. 


1 


1 Introduction 


Gauge theories in two spatial dimensions may contain a Chern-Simons term |[l|— 
[fj] which explicitly breaks parity and time-reversal symmetry. In these theories 
the gauge held acquires a mass and the charged particles obey fractional anyonic 
statistics Non-Abelian Chern-Simons theories have intricate relations to knot 

theory, the Jones polynomials [Toj] , and to 2-dimensional conformal held theories 


11 , 12 


Abelian Chern-Simons theories have been used to facilitate bosonization 

Furthermore, Chern- 


or Fermi-Bose transmutation in (2 + 1) dimensions 13-116 


Simons gauge theories are of central importance in the context of the fractional 
quantum Hall effect 17-19] and other condensed matter systems 20-2^| • They also 


play an important role for topological quantum computation 23 


28]. By braiding 


world-lines of anyonic quasi-particles, one can accumulate appropriate non-Abelian 
Berry phases 2Q, [30], which can encode quantum information. When the anyons 


obey a sufficiently complex version of non-Abelian braid statistics, they can be em¬ 
ployed to realize the quantum gates that are sufficient to realize universal quantum 


computation [31I-I33]. The idea of topological quantum computation is attractive 


because the quantum information is then naturally protected from decoherence by 
the topological nature of the non-Abelian Berry phases. In particular, information 
is not stored locally but is distributed throughout the entire system. The toric code 
is a (2 + l)-d Z(2) lattice gauge theory which can be used as a topologically pro¬ 
tected storage device for quantum information f23j. This theory has charges and 
dual charges with mutually anyonic statistics. 

In this paper, we derive extensions of the toric code from a doubled compact 
lattice Chern-Simons-Maxwell theory with Abelian gauge group 17(1). Similar to 
fermions, topologically massive gauge fields also suffer from a lattice doubling prob¬ 
lem. Here we are not trying to circumvent this problem but work with a lattice 
gauge held and an independent gauge held associated with the dual lattice. The 
fundamental gauge degrees of freedom are then associated with a cross formed by 
a link and its corresponding dual link. The Chern-Simons term couples the two 
lattices and implies that the canonically conjugate momenta (i.e. the electric held 
strengths) of the original and the dual gauge held do not commute. In ordinary 
lattice gauge theories bb the held algebra is link-based. This framework has 
also been used in studies of Chern-Simons gauge theories on the lattice B> 381-42 


In our lattice formulation of a doubled Chern-Simons gauge theories, on the other 
hand, the held algebra is cross-based. Such a system was hrst investigated in 
Here we concentrate on the relation of this theory to the toric code. 
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In ordinary lattice gauge theory with a link-based held algebra every link has a 
“mechanical” analog. It behaves like a “particle” moving in the group space. For 
example, the dynamics of a link variable in a compact Abelian 17(1) lattice gauge 
theory is analogous to the one of a quantum mechanical particle moving on a cir¬ 
cle. Similarly, in our unconventional lattice gauge theory with a cross-based held 
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algebra, the “mechanical” analog of each cross is a charged “particle” moving on a 
2-dimensional group space torus U( l) 2 threaded by an abstract “magnetic” field [ii3 - 


461 ]. The Dirac quantization condition for the abstract “magnetic” flux then implies 
the quantization of the level k — the prefactor of the Chern-Simons term. Inter¬ 
estingly, the corresponding cross-based Hamiltonian has two self-adjoint extension 
parameters, which naturally enter the quantum theory as external parameters, while 
the classical theory is insensitive to these parameters. Remarkably, the self-adjoint 
extension parameters (which are associated with the links and the dual links) them¬ 
selves form two non-dynamical U( 1) lattice gauge fields. This reduces the manifest 
gauge symmetry of the quantum theory from 17(1) to Z(k). This quantum mechan¬ 
ical breaking of the gauge symmetry could be called an “anomaly”. However, it is 
more like the explicit breaking of CP invariance caused by a non-zero 0-vacuum an¬ 
gle in 4-dimensional non-Abclian gauge theories. It is intriguing that the quantized 
doubled Chern-Simons-Maxwell lattice gauge theory dynamically reduces its mani¬ 
fest gauge symmetry from 17(1) to the Z(k) gauge group of an extended toric code. 
The full dynamics reduces to the one of the toric code in the pure Chern-Simons 
limit of infinite “photon” mass. 

The main purpose of this paper is to perform a detailed mathematical derivation 
of an extended version of the Z(k) variant of the toric code as a limit of doubled 
17(1) Chern-Simons-Maxwell lattice gauge theory. We put particular emphasis on 
the role played by the self-adjoint extension parameters of the cross-based Hamilto¬ 
nian. This embeds the toric code in a wider theoretical framework and thus provides 
a broader perspective on quantum information processing. In particular, it would 
be interesting to investigate whether Berry phases that arise from adiabatic changes 
of the external self-adjoint extension parameters can be utilized for quantum in¬ 
formation processing. In this paper, we do not yet address these questions. We 
also view the present paper as a first step towards similar studies in the context of 
non-Abelian Chern-Simons theories on the lattice. Variants of the toric code with a 
discrete non-Abclian gauge group 47], which are discussed in the context of topo¬ 
logical quantum computation, may be related to Chern-Simons gauge theories with 
continuous gauge groups in a similar manner. 


The rest of the paper is organized as follows. In Section 2 we discuss Chern- 
Simons-Maxwell gauge theory of a single Abelian gauge held in the continuum, 
while in Section 3 we investigate the doubled theory with two Abelian gauge fields. 
In Section 4 we regularize this theory on a lattice and its dual lattice, with a cross- 
based held algebra using non-compact Abelian lattice gauge helds. In particular, 
we discuss the mutual anyonic statistics of the charges moving on the original and 
the dual lattice. In Section 5 we turn to compact Abelian gauge helds, which leads 
to the quantization of the level k as well as to the generation of a non-dynamical 
background lattice gauge held formed by the self-adjoint extension parameters of the 
cross-based Hamiltonian. In the limit of a large “photon” mass this theory reduces 
to the toric code. Finally, Section 6 contains our conclusions. 
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2 Chern-Simons-Maxwell Theory of a Single 
Abelian Gauge Field in the Continuum 


In this section we investigate Abelian Chern-Simons-Maxwell theory in the contin¬ 
uum. After considering a single Abelian gauge held, in the next section we will 
study a doubled theory which will turn out to arise in the continuum limit of the 
lattice theory that we discuss later. 


2.1 Lagrangian and Hamiltonian 


Let us consider a (2+l)-d Abelian gauge field A^x) with the Chern-Simons-Maxwell 
Lagrangian 

+ ( 2 . 1 ) 

where e is the electric charge and the metric is = diag(l, — 1, — 1). The field 
strength is given by F liu (x) = d^A^x) — d v A^(x). Without loss of generality, in the 
following we assume that the prefactor of the Chern-Simons density k > 0. It is 
interesting to note that the Chern-Simons density is not gauge invariant. However, 
its variation under a gauge transformation 


A'(x) = Ap{x) - dpxix), 


is a total divergence 


sc = a -c = -&• 


4)r W 


XX.V 


( 2 , 2 ) 


(2.3) 


It should be noted that the Chern-Simons term explicitly breaks both time-reversal 
and parity (i.e. the reflection on a spatial axis). In order to derive the corresponding 
Hamiltonian, we now fix to the temporal gauge A 0 (x) = 0, which implies 

C = -^rAf - ^-UiAAA 2 + —etjAiAj. 

2e 2 1 2e 2 v J An 3 3 

The momentum canonically conjugate to Ai is then given by 

5C 

~ ~ b2 J 




1 • k 

= —Ai(x) + —e ij A j (x), 


5Ai(x) e 

such that the classical Hamilton density takes the form 

\2 


"H — TliAi C— 2e2 A + 2 e 2 

2 


■l 1 H 
2 


n,; 


k_ 


1 

2^ 2 ' 


+ ir^i^ijdiAj) — —Ef + 


e 

~2 


2e 2 


Here we have identified the electric and magnetic fields as 

k 


Ei(x) = n*(x) - —e ij A j (x), B{x) = diA 2 (x) - d 2 A 1 {x). 

47r 


(2.4) 

(2.5) 


( 2 . 6 ) 


(2.7) 
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2.2 Solutions of the Classical Equations of Motion 


Before we quantize the theory, we consider its classical equations of motion 

he 2 

d^F^(x) + —e vpa F^{x) = 0, (2.8) 

47T 

which in components take the form 

1 ke 2 k 

Ei(x ) + —(;it), F>(.r) + —e ij E j (x) = 0, diE^x) + — B(x) = 0. (2.9) 

e- An An 

The second equation is the Gauss law. In particular, the magnetic held acts like a 
charge density. In addition, the held strength obeys the Bianchi identity 

e pvp d^F vp {x) = 0. (2.10) 


Let us make the simple plane wave ansatz 

Ei(x) = Ci cos (p- x — ut) + Di sin(p • x — ut). (2.11) 

Inserting this in the equations of motion, one obtains 

. _ ^ e 2 m m 

u = y/M 2 + p 2 , M = —, Ci = ceijPj - d — Pi, Di = deijPj + c — p h (2.12) 

An lo co 

where c,dGl are arbitrary constants. We have identified M as the topologically 
generated “photon” mass of the Abelian gauge held, which is proportional to the 
prefactor k of the Chern-Simons term. In the absence of the Maxwell term, i.e. 
when e 2 —» oo, the gauge held becomes infinitely heavy. In the pure Maxwell theory 
(with k — 0), on the other hand, the gauge held remains massless, and plane waves 
are purely transverse. 


In order to quantize the theory, we now impose canonical commutation relations 

[n* (x ), Aj(y)\ = -i6ij6(x - y ), (2.13) 

which imply the following commutation relations between the electric and magnetic 
hclds 

k 

[■ E i(x ), Ej(y)\ = -i—e i:j 5(x - y), 

[B{x),B{y)] = 0, 

[Ei(x),B(y)\ = ie ij d j S(x - y). (2.14) 

The derivative in the third equation is with respect to x (not y). As a consequence 
of the Chern-Simons term, the two components of the electric held do not commute 
with each other. It is easy to check that the Hamiltonian 

h = + P' 15 ) 
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commutes with the infinitesimal generators of local gauge transformations 

G(x) = diEiix) + [G(x), H] = 0. (2.16) 

Z7T 

As usual in a gauge theory, physical states | T) must obey the Gauss law G(x) | \&) = 0. 
Upon quantization, the classical plane wave solutions then turn into free “photon” 
states of mass M = tF. 


3 Doubled Continuum Theory with two Abelian 
Gauge Fields 


In the next section, we will regularize Chern-Simons-Maxwell theory on a particular 
lattice, which will result in a doubling problem, similar to the well-known lattice 
fermion doubling problem. The continuum limit of that lattice theory is a dou¬ 
bled Chern-Simons-Maxwell theory with two Abelian gauge fields, which we first 
investigate directly in the continuum. 


3.1 Lagrangian and Hamiltonian 


We now consider two Abelian gauge fields A /1 (x) and A^x) with the Lagrangian 

£ = ~^ F ^ F " U ~ 4p^ {A^Ap + ArdTAP). (3.1) 

Without loss of generality one can assume that the dual charge equals the original 
charge (i.e. e — e). If this is not the case a priori, one can achieve this by a simple 
rescaling of the fields. Unlike the theory with a single gauge field, the doubled 
theory is invariant under time-reversal and parity, because one can treat A^ix) as a 
pseudo-vector. 

Again, fixing to the temporal gauge Aq{x) = 0, A Q {x) = 0, one obtains 


C - —A] - —{eijdiAj ) 2 + —A i - —( e^diA^Y + — e ij (A i A j + AiAj), (3.2) 


which yields the following canonically conjugate momenta 



5C 

1 • . . 

k ~ 

Ilj(x) 

sAi(x) 

= -zAi(x) 
e z 

+ ^ e v A j( x - 


5C 

1 ~ 

k 

n. t (x) 

SAAx) 

= —A(x) 

+ y:''o 


(3.3) 
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The classical Hamilton density then takes the form 


TL — TliAi + II jHj — C — 


A 2 + -^{eijdiAj ) 2 + + 


2e 2 


2e 2 


= - n 


k_ j 
4Tr eijAj 


+ ^(eydiA) 2 + y 


n, 


A 

47 r e *i A i 


2e 2 

2 


tijdiAj ) 


+ 2 


2 Ei + 2e 25 + 2 + 2e 2 


H 2 . 


(3.4) 


In this case, the electric and magnetic fields are given by 


Ei{x) 

= AA) 

k ~ . . 

- 4 -eijAjix), 

B{x) 

= d 1 A 2 {x) 

- d 2 A l (x), 




k 





Ei{x) 

= Ilj(x) 

~ 4 -^j A j{ x h 

B(x) 

= di A 2 (x) 

- d 2 A 1 (x). 

(3.5) 


3.2 Solutions of the Classical Equations of Motion 


The classical equations of motion of the doubled theory are 


d^F^x) + —e vp<J F^{x) = 0, d^F^x) + —e vpa F^{x) = 0, (3.6) 

47r 47T 

which in components take the form 

1 ke 2 ~ k ~ 

Ei(x ) + —e ij d j B(x) + —e ij E j (x) = 0, d\Ei(x) + —B(x) = 0, 

e 2 7T 27T 

~ 1 ~ ke 2 ~ k 

Ei(x ) + —e ij d j B(x) + —eijE^x) = 0, diE^x) + — -B(z) = 0. (3.7) 

As before, we make the plane wave ansatz 


E'j(x) = C cos(p • x — Lot) + Di sin(p • x — cot), 
Ei(x) = Ci cos(p • x — cot) + Di sin(p • x — Lot), 


which again implies lo 


A/M 2 + p 2 , and M — as well as 


Ci = 



< < i.jp.j - d — pi, 

A 


LO 



M 


Ci = 

C^ijPj d Pi-) 

CO 

Di 


<l< ijPj + c— Pi, 

CO 

7 M 

deijPj + c—pj. 


LO 


(3.8) 


(3.9) 


In this case, there are four independent amplitudes c, d,c,d E R, and thus there are 
two independent “photon” modes, both with the same mass M. 
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3.3 Canonical Quantization 


Canonical quantization now amounts to 

[n 4 {x),Aj(y)] = [H i(x),Aj(y)\ = - 18 ^ 8 (x - y ), 

[n,; (x ), Aj{y)\ = [fii(x), Aj(y)\ = 0, (3.10) 

which imply the following commutation relations between the electric and magnetic 
fields 

[Ei(x),Ej(y)] = [Ei(x), Ej(y)\ = -i^- €ij S{x - y), 

[Ei(x), Ej(y)] = [Ei(x), Ej(y)} = 0, 

[B{x),B{y)\ = [B{x),B{y)\ = [B{x),B{y)\ = 0, 

[Ei(x), B(y )] = [Ei(x), B(y )] = ie^d^x - y), 

[E i (x),B(y)] = [E i (x),B(y)] = 0. (3.11) 

The Hamiltonian now takes the form 

H = / ^ (j e - + ^ + P- + w s2 ) ■ (3 ' 12) 

It commutes with the infinitesimal generators of two local gauge transformations 

G(x) = diEi(x)+-^-B(x), G(x) = d\Ei(x)+-^B(x), [G(x),H] = [G(x),H] = 0. 
Ztt Ztt 

(3.13) 

The Gauss law now constrains physical states to G(x)| v h) = G(x)| v h) = 0. The 
classical plane wave solutions then turn into the states of two free “photons”, both 
of mass M = 

Z7T 


4 Non-compact Chern-Simons-Maxwell Theory 
on the Lattice 


In this section, we regularize Chern-Simons-Maxwell theory on a particular spatial 
lattice. We work in the Hamiltonian formulation and thus leave time continuous. 
As we will see, the lattice theory suffers from a doubling problem and thus reduces 
to the doubled Chern-Simons-Maxwell theory in the continuum limit. 


4.1 Cross-based Degrees of Freedom 

For simplicity, we consider a square lattice. The generalization to other lattice 
geometries is straightforward but not very illuminating. Most previous lattice con¬ 
structions of Chern-Simons theories have placed all dynamical gauge degrees of 




Figure 1: [Color online] a) Square lattice A (squares, solid links) together with its 
dual lattice A (diamonds, dashed links). The gauge fields A Xji and A x>i are associated 
with two links forming a cross centered at the point x. The cross centers x do not 
belong to the lattices A or A, but form two other sublattices X (filled circles) and X 
(open circles), b) Two types of crosses of link-pairs, A Xjl ,A X:2 and A Xj i,A X j2 , whose 
cross centers x belong to the two distinct sublattices X and X, respectively. 

freedom on the same lattice. Here, similar to the lattice geometry used in [43|, we 
simultaneously place independent gauge degrees of freedom on the original as well as 
on the dual lattice. This is natural from a geometric point of view, but it inevitably 
leads to the doubling of the massive “photon” mode. This is no problem in the 
present context, because it is our goal to explicitly connect lattice Chern-Simons- 
Maxwell theory to the toric code, which indeed requires to start from the doubled 
theory. 

In standard lattice gauge theory, the basic gauge field algebra is link-based, 
i.e. the gauge field variables as well as their canonically conjugate electric field 
operators are defined on a link. In particular, field operators residing on different 
links commute with one another. As we have seen, in the doubled Chern-Simons- 
Maxwell theory the canonically conjugate variable to A; is n, ; = -fA,- + Xe^Aj. For 
the corresponding lattice theory, this implies that the gauge variables associated 
with different links no longer necessarily commute. As we will see, it is natural to 
associate the non-commuting gauge field variables with a link and the corresponding 
dual link. Consequently, the basic gauge algebra of doubled lattice Chern-Simons 
theory is not based on a single link, but on a pair of links which are dual to each 
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other. The original square lattice together with its dual lattice is illustrated in Fig 03 

Let us consider a link and the corresponding dual link which form a cross centered 
at a point x. Note that x is neither a site of the original nor of the dual lattice, 
but just an intersection point of the link and its dual link on a sublattice A" or 
X (cf. FigOJ- We denote the non-compact vector potential on the original link by 
A Xj i. The vector potential on the corresponding dual link is then given by 
It should be noted that for a given cross, the direction i of the original link variable 
A Xji and the orthogonal direction j of the dual link variable A x j are determined by 
the location of the intersection point x on sublattice A" or X. 


4.2 Lagrangian and Hamiltonian 


In the continuum the integrated Lagrange density (i.e. the Lagrange function) is 
given by 


L= d~x 


1 i2 1 D 2 


1 A 2 


k 


2e 2 A ‘ - ^ B ~ + A ‘ - ^ B ~ + ^‘Mi A l + A i A i) 


Correspondingly, on the lattice we obtain 


(4.1) 


L = J2 a2 

xex 


1 12 IV k 


E 

xex 


a 


2 e 2 

1 


Ar,i + -^A x2 + —(A x> 1 A Xj2 - A X ^A X ^) 


47T 


Id 2 k A 


Az,2 + 7^2 ^x,l + ^jl(4r,l4r,2 - A X ^A X ^) 


2 e 2 a ' 2 2e 2 ' 


- E^-E“ 2 A« 


cGA 


ce A 


2 e 2 


(4.2) 


In the sums over the sublattices A", A", A, A the factors a 2 (where a is the spacing 
of the original or the dual lattice) arises as the area of an elementary plaquette. On 
the lattice, the magnetic fields are defined as 


B x 

B x 


a i^-i 2 ’ 1 


A 


■*—§ 2,1 


+ Ac+§1,2 Az+§2,1 Ae—§1,2 J ) 

+ Ae+§1,2 — Ae+§ 2,1 — Ar-fi,2 J ) 


(4.3) 


with i being the unit-vector in the f-direction. Note that the magnetic field B x , 
whichms built around a plaquette of the original lattice, is indexed by the dual site 
x G A at the center of this plaquette. The magnetic fields are invariant against 
lattice gauge transformations 






(4.4) 
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On the lattice, the canonically conjugate momenta are given by 


n x ,i 




dL 2 { 1 ; k 7 

- — a I ~ A x i + — 6 ijA. 

dA x>i \e 2 4tt 
dL 2 ( 1 T k . 

^ 61 I o -^-xA “1“ , 

dA x/l V e 


2^ x ’ 1 T 4 7 r c o y ^o 


(4.5) 


Note that on the lattice Ti x ^ has a different dimension than II,; (re) in the continuum. 
The corresponding classical Hamilton function then takes the form 

h = + na,,2^4z, 2 ) + y ] (n x , 2A Xi 2 + n x ,iA Xi i) — l 


xex 

e _ 2 
2 


x&X 


Y ° 2 ( E l,i + E l?) + 77 'Y + E h) 


xGX 


1 E^+^E^: 


xex 

2 


2 e 2 


x£A x£A 

In this case, the electric fields are given by 

1 k ~ — 1 — k 

a 2 47 t a 2 47 r 


(4.6) 


(4.7) 


4.3 Solutions of the Classical Equations of Motion 


Before we quantize the theory, let us first consider the classical equations of motion 
on the lattice, which in components take the form 

Y “ E x-fi,i) + l^ E x = 0 , 
i 

1 1 ke 2 — 

E x,i + —^€ij — (B x+ aj — B x — aj ) + ——6 lJ E x j = 0 , 


1 


Eke 


a ' x+ * i ’ i x ~^ 


2 J ' 27r 

E x -<i.i) + 7~B x = 0 , 


2 tt 


~ 11 — — ke 2 

Ex ’ i E ^ 2 e v q (-^®+§ 5 ' — E x-%j) A ~^ e ij E x,j = 0 - 


(4.8) 


The first equation corresponds to the lattice Gauss law, which is illustrated in FigJ2f,. 
The geometry underlying the second equation is depicted in Figj2f>. As in the 
continuum, we make the plane wave ansatz 


E Xy i = Ci cos (p- x — cot) + D, sin(p • x — cot), 
E X j = Ci cos {p ■ x — cot) + D, sin(p • x — cot). 


(4.9) 
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Figure 2: [Color online] Geometry of the lattice equations of motion: a) The Gauss 
law arises from the divergence of the electric field, i. e. from the electric fluxes flowing 
out of a lattice point x, and from the magnetic field associated with a plaquette formed 
with the dual links, b) The other two equations of motion relate the time-derivative 
of an electric field E x i on a link centered at x to the dual electric field (i 3 E XJ as well 
as to the difference of magnetic fields residing on two nearest-neighbor plaquettes. 


Inserting this in the lattice equations of motion, we obtain the dispersion relation 
co = sjM 2 + p 2 , M — Pi = - sin^p. (4.10) 

Z7T a 2 j 

Note that the mass of the gauge field is the same as in the continuum, while the 
energy-momentum dispersion relation is distorted by lattice artifacts for lattice mo¬ 
menta p close to the edge of the Brillouin zone. In analogy to the continuum result, 
we obtain 



A 

, A AlA 

C^ijPj d Pj, 

CO 

= denPi + c — p, 

CO 

~ A M 

___ 

T A M 

ceijPj - d — Pi, 

CO 

A 

= deijPj + c— p, 

UJ 


(4.11) 


As in the continuum, there are four independent amplitudes c, d,c,d G R, and thus 
the massive “photon” mode is again doubled. 


4.4 Canonical Quantization 


Upon canonical quantization we now obtain 




The resulting commutation relations for the vector potential and the electric field 
are 


[A X ,il Ay t j\ 

\E X ,il Ayj\ 

[E X ,il Ayj\ 

[E X ,i, Ey,j] 
[E x ,ii Ey,j] 


[A X ,il Ayj] 

[E X ,i, A yj] 
[E X ,i, Ayj 
[E X ,il EyJ 
[E X ,il EyJ 


[^-x,n 0, 

1 


Auij ^ °xy, 


] = 0, 


■— J_A 

'W xy ' 


= 0. 


(4.13) 


In the continuum limit a —> 0, jp8 xy approaches 8{x — y). It is interesting to note 
that the commutation relations between A x j and e l3 A XJ with E xy and e l3 E x j, which 
reside on the cross formed by a link and its dual link, correspond to those of a charged 
particle moving in a 2-d plane in a constant magnetic field E, j n this “mechanical” 
analog, A x i and e.jA x j play the role of the particle’s x- and ^-coordinate, while E x l 
and 6ijE x j represent the corresponding momenta. Like for a charged particle in a 
magnetic field, the commutator of the momenta is non-zero and given by the abstract 
“magnetic” field d_. In the next section, we will compactify the gauge field from R 
to 17(1), which will naturally lead to a quantization of the abstract “magnetic” field 
in integer units k G Z. For the moment, however, k is not quantized. 


The commutation relations between electric and magnetic fields now take the 
form 

[E x j, By] = [E x j, By] = ieij—(5 + aj — 8 a A), 

(1 A A 

[E X J, By] = [E x>i , By] = 0, 

[B x ,By] = [B x ,By] = [B x ,By] = 0. (4.14) 


The Hamiltonian of eq. 04.61) commutes with the infinitesimal generators of local 
gauge transformations, both on the original and on the dual lattice 


G x _ 1 ( 

a 2 ^ E x+ 


E x -±lj) + ^ X E 


G 

a* 




+ [G„J/] = 0, xeA. (4.15) 


Note that G x has a different dimension than G(x) in the continuum. In general, 
generators of gauge transformations at different sites of the same sublattice com¬ 
mute with each other, i.e. [G x , G y ] = 0, [G x , G y ] = 0. Let us also compute the 
commutation relations of the generator of a gauge transformation on the original 
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Figure 3: [Color online] Consistency of the Gauss laws on the original and on the 
dual lattice. The gauge generators G x and G x (with x = x + |1 + |2j commute, 
i.e. [G X ,G X ] = 0. This is non-trivial, because, e.g., the electric fields E x+ and 

E x _a ,2 2 that contribute to G x and G x do not commute since they belong to the same 
cross formed by a link and its dual link. 


lattice and a neighboring site on the dual lattice (cf. FigJS]) 

[G x ,G x+ al + af\ = —a [E x+ al 1 , E x+ al^ 2 ] — a [E x+ a 2 t 2 , E x+ |2,l] 

+ ~2ff^ x + §i,i’ i+fd + ~^ff[Ex+%2,2 > -^x+§i+§2] 

+ ^[4+52,1. G + ^A+fi, 2, S,]=0. (4.16) 

Hence, generators of gauge transformations at different sites always commute. 


The vacuum state |0) is gauge invariant and satisfies the Gauss law G x |0) = 
(j x |0) = 0. Upon quantization, the classical plane wave states turn into the states 
of two “photons”, each of mass M = (f-, with the lattice dispersion relation E(p) = 
\JM 2 + p 2 , Pi — - sin 2^. A state | Q, Q) that contains static charges Q x E Z (with 
x E A) and static dual charges Q x E Z (with x E A) satisfies 

G X \Q,Q) = Qx\QiQ), x E A, G X \Q, Q) = Q X \Q, Q), x E A. (4-17) 


4.5 Mutual Statistics of Original and Dual Charges 

Let us consider Wilson’s parallel transporter U Xii = exp^aA^), which moves a 
charge from x — |z to x + |b Linder the unitary transformation 

V = JJ 14 = JJ exp(iXxG x ), (4.18) 

xGA xGA 
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Figure 4: [Color online] A Wilson loop Wc on the original lattice encircles a set of 
dual charges Q x (filled diamonds). Up to exponentially suppressed corrections due to 
the finite “photon” mass, the vacuum expectation value of the Wilson loop depends 
only on the total encircled dual charge, not on further details of the closed loop C. 


which implements a general gauge transformation in Hilbert space, the parallel trans¬ 
porter transforms as 


K,i = VU x ,iV ] = exp(ix x _a-)U Xii exp(-ix x+ ai), xeX,X. (4.19) 

The corresponding commutation relation takes the form 

[Gy, U x ,i] = - hi - n y-fhi’ ^ 4 ] = (fiy+fi,x - U x ,i. (4.20) 


Let us now act with U x ^ on the charged state \Q,Q), with external charges Q x at 
x G A and Q x at x G A (cf. eq. fl4.17p ). in order to investigate its effect on the charges 


G y \Q',Q) = G y U x>i \Q,Q)^U Xii G y \Q,Q)+(S y+ a lx -S y a lx )U x>i \Q,Q) 


- (Qy + 8 y+ a- x ~ S y _a~ ijX ) \Q', Q) 


(4.21) 


As expected U X)i acts as a charge transport operator which moves a charge from 
x — |i to x + | i, i.e. the state | Q', Q) = U X ^\Q, Q) has charges 

Qy = Qy + fiy+fi,x ~~ (4.22) 


Let us now transport a charge around the closed loop C (on the original lattice), 
so that it returns to its initial position (cf. Figj4]). The wave function then turns 
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into 


(4.23) 


W C \Q,Q)= I] U Xji \Q,Q) = exp(i$ c )\Q,Q), 

( x,i)&C 

i.e., it picks up a phase 

= £ “ 2 B, = V yG, - £ (4.24) 

xgScA x€S cA c 

which is given by the magnetic flux through the surface S that is bounded by C — dS. 
Here x denotes a point on the dual lattice at the center of a plaquette on the original 
lattice which belongs to the surface S. Eq. fl4.24p results from Stoke’s theorem on the 
lattice as well as from eq. (l4.15p . The first term on the right-hand side of eq. (I4.24H 
counts the dual charges Qs = J2 x gsca ^ x enc i r cled by the loop C. In a pure Chern- 
Sirnons theory (without a Maxwell term) this would be the only contribution. In 
that case, a charge that is transported around a closed loop C on the original lattice 
picks up a topological phase ex.p(2mQ$/k). In other words, charges on the original 
and on the dual lattice have mutual anyonic statistics, with the statistics angle 
The topological phase is an adiabatic Berry phase. However, in contrast to generic 
Berry phases, the topological phase <f>e is even insensitive to the shape of the curve 
C, as long as it encircles the same total charge Qs on the dual lattice. Obviously, 
corresponding rules about charge transport also apply to the dual lattice, with U x/l 
being replaced by U X)l = exp {iaA x Q. In the Chern-Simons-Maxwell theory that we 
investigate here, a second term arises on the right-hand side of eq. (l4.24p . This term 
is due to the “photon” cloud surrounding the dual charges. Since the “photon” is 
massive, the effect of this term is exponentially suppressed at large distances 0. 
Still, it slightly affects the topological nature of the Berry phase and causes a certain 
sensitivity to the precise location of the curve C. The central result of mutual anyonic 
statistics of original and dual charges with statistics angle ^ remains valid, as long 
as the charges are moved around each other at distances much larger than the size 
of their massive “photon” clouds. 


5 Compact Chern-Simons-Maxwell Theory 
on the Lattice 


In the previous section, we have considered a non-compact lattice theory with gauge 
group R. As we have seen in eq. fl4.13p . the local degrees of freedom on a cross formed 
by a link and its dual link then have a “mechanical” analog, namely a charged 
“particle” moving in the 2-d plane R 2 in the background of an abstract “magnetic” 
field A. Now we will compactify the gauge field and thus work with the gauge group 
1/(1). The “mechanical” analog then corresponds to motion restricted to the torus 
U( l) 2 . This leads to the quantization of the level k in integer units, i.e. k e 7L. As we 
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will see, the electric field part of the Hamiltonian has a 2-parameter family of self- 
adjoint extensions. Interestingly, fixing the two self-adjoint extension parameters 
explicitly breaks the U( 1) gauge symmetry down to Z (k). 


5.1 Commutation Relations of the Compact Theory 


In the compact theory, A x and A X j turn into phases of the parallel transporters 


U Xji = exp (iaA X}i ) = exp (ip x ,i) G 17(1), U X)i = exp (iaA X:i ) = exp (ip X j) G 1/(1), 

(5.1) 


which now are the truly fundamental degrees of freedom, and 

k ~ ~ k 

(^E x j i ~a Q*E x i i ~t ( 5 - 2 ) 

’ 47T ’ 47T 

Based on eq. 04.13p . the local commutation relations now take the form 


[U X ,l, Uyj] 
[E X ,i, Uyj] 
[E X j, Uyj] 
[E X j, Eyj] 
[E X j, Eyj] 


[U X j, Uyj] [U x j, Uyj] 0, 

&ij&xyUx,ii \E x j, Uyj] $ijS X yUx,ii 

CL CL 

\E X j, Uyj] 0, 

~ k 

[E x ,i, Eyj] = ~^2 na 2 

\E X ,il Eyj] = 0. 


(5.3) 


Since, unlike A X j and A x>i , p x , t and p X j are angles, they do not represent self-adjoint 
quantum mechanical operators. 


5.2 Self-adjoint Extensions of the Electric Contribution to 
the Hamiltonian of a Single Cross 

Let us consider a single cross, formed by a link and its dual link. For concreteness, 
we consider a cross x G X , i.e. the link variable p is associated with the 1-direction, 
while p is associated with the 2-direction. The situation for x G X is completely 
analogous. In order to keep the notation simple, in this subsection we suppress the 
link indices (x,i), which are uniquely determined for a specific cross x G A". The 
“mechanical” analog of the problem then corresponds to a charged “particle” with 
“spatial” coordinates (p, p) moving in the 2-dimensional group space subject to the 
cross-based Hamiltonian 

H + = + El 2 ) = ( ^fu-A + E 2 ). (5.4) 
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It should be noted that, after compactification, a dual charge e can no longer be 
turned into e by a field redefinition. In order to maintain the symmetry between 
the original and the dual lattice, we simply put e — e. The electric field operators 
are now given by 


aE = — id^ + a((p, 0, aE = + a((p, 99 ). 


(5.5) 


Here (a, a) is an abstract vector potential on the group space torus U( l ) 2 (parame¬ 
trized by (<p,0)) with 

a(<P,<P) = —^<P, = (5.6) 

which gives rise to the abstract “magnetic” field 

k 

%, <p) = d v a(ip, 0) - d$a(<p, 0) = — ■ (5.7) 

The level k thus determines the value of the constant abstract “magnetic” field 
b((p,0) (which should not be confused with the actual magnetic fields B x and B x ). 


Since the electric field part of the Hamiltonian E[ + contains E 2 and E 2 , these 
operators must be self-adjoint. For mathematical details related to the theory of 
self-adjoint extensions we refer to 481 4911 . The problem of finding the self-adjoint 
extensions of H + has been discussed in [46] for the “mechanical” analog of a charged 
particle moving on a torus in a constant magnetic field. The same mathematical 
solution applies here. The self-adjoint extension parameters enter the boundary 
conditions for the wave function ^( 99 ,^) on the group space torus U( l) 2 . This 
wave function will extend to the wave functional of the lattice field theory, once 
we combine all crosses to form the entire lattice. For the moment, we continue to 
consider a single cross in isolation. 


While the abstract “magnetic” field b is constant, and thus obviously periodic 
over the group space torus, the abstract vector potential (a, a) obeys the boundary 
conditions 

a(ip + 2ir, 0) 
a(</?, 0 + 27r) 
a((p + 27 t , 0 
a((/?, </9 + 27 t) 

Let us introduce the transition functions 

(*(&) = (5-9) 


k ~ k k 

= ” 4^-2 = “ 2 ’ 

k k ~ k 

= jy + 2 = + 2. 

k 


(5.8) 
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where 9,9 G [0, 27t[ will turn out to be two self-adjoint extension parameters. We 
can then express the boundary conditions of eg. fl5.8[) as twisted periodic boundary 
conditions over the group space torus U( l) 2 

a(tp + 2n, tp) = a(ip, tp) - d^a(tp), 
a{<p, tp + 27 t ) = a(<p, ip) - d^ip), 
a(tp + 27 t , tp) = d(tp, ip) - d$a(tp), 

a{tp,tp + 27 t ) — a(tp,tp) - dpa(ip). ( 5 . 10 ) 

As we see, a and a play the role of abstract gauge transformations (not to be con¬ 
fused with the original lattice gauge symmetry). In order to respect the abstract 
gauge structure, the wave function must obey consistent twisted periodic bound¬ 
ary conditions, i.e. 

T(<^ + 27t, tp) = exp(za(<p))T(<p, tp) = exp ^~^P + ^ V)i 
^(ip,tp + 2n) = exp(i5(<p))d'((p, tp) = exp + i9^j \1 '/(ip, tp). (5.11) 


Note that this equation holds for a cross x G X, while for x G X the 9 and 9 terms 
change sign. Now we see explicitly, that 9 and 9 enter the boundary condition of 
the wave function, which is typical for self-adjoint extension parameters. The above 
relations immediately imply 

T(<^ + 27 t, tp + 27 t) = exp -i— (<p + 27 t) + id j T(<p, <p + 27 t) 


= exp 


k ~ ~ k 

-i-{<p + 27t) + i9 + i-tp + i9 ) \l/(</?, tp), 

k 


^>(tp + 27t, tp + 27t) = exp {tp + 27t) + i9 J T(<p + 27 t, tp) 

/ k k 

= exp ( i-(ip + 27t) + i9 - i-tp + i9 j ^(tp, ip). (5.12) 

Consistency of the two expressions requires exp(27rifc) = 1, which leads to the quan¬ 
tization of the level k G Z. This quantization condition implies that the abstract 
“magnetic” flux through the group space torus U( l) 2 is quantized as well 


/‘27T /»27T 

/ dtp dtp b{tp, tp) = 2nk. 

'0 Jo 


(5.13) 


The quantization of the level k thus manifests itself as a Dirac quantization condition 
of the abstract “magnetic” flux that threads the group space torus. 

Besides the field strength b , the abstract gauge held (a, a) on the group space 
torus U( l) 2 has two additional gauge invariant quantities — the Polyakov loops that 
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arise due to the non-trivia! holonomies of the torus 

/*27r ^ 

(f>(p) = dp a((p, ip) + a (ip) = -kip + 0, 

Jo 

~ r 27r 

4>{ip) = / dp a(p, p) + a(p) = kp + 6. (5-14) 

Jo 

Unlike the classical theory, the quantum theory is sensitive to the complex Aharonov- 
Bohm phases defined by the Polyakov loops 

exp [i(j>{!p)) = exp (—ikp + id), exp (i(f>(p)) = exp (ikp + id). (5.15) 

Remarkably, the Polyakov loops explicitly break the continuous U( l) 2 translation 
invariance of the torus down to the discrete translation group Z(A;) 2 , since the cor¬ 
responding Aharonov-Bohm phases are invariant only under shifts of p and p by 
integer multiples of G-. This corresponds to an explicit breaking of the original and 
dual lattice U( 1) gauge symmetries down to Z (k). 


5.3 Fate of the Original and Dual U( 1) Gauge Symmetries 

As we have just seen, when we fix the values of the self-adjoint extension parameters 
6 and 6 , the original and dual continuous U( 1) gauge symmetries of the classical 
compact doubled Chern-Simons-Maxwell theory are explicitly broken down to dis¬ 
crete Z(k) gauge symmetries. Since gauge symmetries just reflect redundancies in 
our theoretical description and are not actual physical symmetries, one may wonder 
what happened at the quantum level to the redundancy associated with the continu¬ 
ous U( 1) gauge symmetries of the classical theory. In order to answer this question, 
let us perform original and dual continuous U( 1) gauge transformations on the link 
variables 


K,i = ex P(^,J =exp(ix x _«i)£4,iexp(-ix s+ «i) 

= exp(i(^ X)i - Xz+f 1 + 

K,i = exp(i^j) = exp(ix x _ai)U x ^exp(—ix x+ ai) 

= exp (i(p x>i - x x+ a- + Xx-fi))- ( 5 - 16 ) 

While the magnetic part of the Hamiltonian is by construction manifestly gauge 
invariant, the electric part is at least not obviously gauge invariant, because E and 
E depend on p and p. Still, in the non-compact theory, E and E turn out to be 
gauge invariant, i.e. they commute with local gauge transformations G and G. Let 
us return to the compact theory and consider a general finite gauge transformation 
that is represented in Hilbert space by a unitary transformation 

V = n exp (iXxG x ) JJ exp(iXxG x ). (5.17) 

xEA 
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It is straightforward to convince oneself that indeed 


vu Xti v ] = U' Xii , VU X ^ = U' Xii , (5.18) 

with U' xi and U' xi given by eq. (l5.16l) . Furthermore, at least at a formal level, 

VE Xii U f = E Xii , VE X ^ = E Xii , VHV f = H, (5.19) 


which seems to suggest that the compact theory is still U( 1) gauge invariant. How¬ 
ever, this is not the case, because V does not leave the domain of the Hamiltonian in¬ 
variant. In other words, the gauge transformed wave functional ^'\<p, <p\ — V^\<p, tp\ 
no longer satisfies the self-adjoint extension condition eq. (l5.1ip . Unlike the wave 
function (which only depends on the link variables ip and p> of a single 

cross), the wave functional ^[(p,(p\ depends on the entire lattice gauge fields [(p,<p\- 
Under a general gauge transformation, one obtains 

'J'V <p\ = VV[(p, tp\ = exp JJ exp ^[(p 1 , &], (5.20) 

xeA ^ n ' xe X V 77 J 


where p>' and ip' are the gauge transformed fields of eq. fl5.lbp . and the plaquette 
magnetic fields are given by 


exp(fa 2 H r ) 

exjp(ia 2 B x ) 


exp(*(¥>*-§2,l + Vx+fl ,2 L Px+ §2,1 VTc-flg))’ x E A, 

exp(i(^_ f 2 ,r + £*+f i, 2 - ^+| 2 ,i - <P X - fig))’ x G A.(5.21) 


When one applies the unitary transformation to the wave function of a single cross, 
one finds that the transformed wave function obeys the boundary condition 


^\<p + 27t, <p) = exp i^cp + ^'((p, <p), 

fl y(<p, <p + 27t) = exp U i^p + iff') ^'(<p, <p), (5.22) 


where the self-adjoint extension parameters 6 and 9, which can be associated with 
the link and its dual link, respectively, have been transformed to 

exp (iff x>i ) = exp (i(0 x>i + kx x+ p - kx x _p)), 

ex P (K,i) = exp (i(0 Xii + kx x+ a -- kx x _a-)). (5.23) 

This means that the self-adjoint extension parameters exp (i9 X j) and exp (i6 Xt i) them¬ 
selves also transform as 17(1) gauge fields. However, unlike U x/l and U x> i, they trans¬ 
port k units of charge. As a consequence, exp (i9 X)i ) and exp (iff X j) are invariant 
against Z (k) gauge transformations. Since a U( 1) gauge transformation changes 
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the self-adjoint extension parameters 9 and 9, and thus the domain of the Hamil¬ 
tonian, it can no longer be considered as a symmetry of the theory. However, the 
redundancy associated with the 17(1) gauge symmetries of the classical theory is still 
present at the quantum level. Continuous 17(1) gauge transformations simply lead 
from one domain of H (characterized by 6 and 9) to a unitarily equivalent domain 
characterized by 9' and 9'. The gauge invariant physical content of the theory is 
thus not determined by the link parameters 9 X and 9 x i themselves, but by their 
plaquette field strengths 

exp(i^) = exp(i(0 ;r _|2 )1 + 9 X+f i 2 - 9 X+ - 9 X _ a ij2 )), x G A, 

exp(i^ s ) = exp(i(0 x _« 2)1 + 0 x+ « i)2i ;v 0 x+ .2q-^-fiy)), x e A. (5.24) 

A natural choice is r/ x — rj x — 0 for all x. Another choice is r/ x = rj x = n. Both of 
these choices leave the charge conjugation symmetry C intact, while other choices 
break C explicitly. The parity symmetry P is already explicitly broken by the 
Chern-Simons term, ft should be pointed out that, unlike U X j and U x ^, exp (i9 X j) 
and exp [i9 X)i ) are not dynamical gauge degrees of freedom, but parameters that 
define a super-selection sector. In this sense, they are analogous to the 9 vacuum 
angle in 4-dimensional non-Abelian gauge theories. In particular, different values of 
9 define different vacuum sectors. For 0 / 0,7r, in non-Abelian gauge theories the 
vacuum angle explicitly breaks the CP symmetry. 

What is the significance of the reduction of the manifest gauge symmetry from 
17(1) to Z {k) at the quantum level? In particular, should one interpret this explicit 
quantum mechanical symmetry breaking as an anomaly? It is important to note 
that the quantum theory is sensitive to external parameters (the gauge fields 9 and 
9 representing Aharonov-Bohm phases) which do not affect the dynamics at the 
classical level. Remarkably, in this case the external parameters are Z (k) gauge 
invariant but change under I/( 1) gauge transformations. 


5.4 Local Magnetic Translation Group on a Single Cross 


As we have stressed before, the actual manifest local symmetry of the Hamiltonian 
in a given domain is reduced to Z(/c). To investigate this symmetry, first let us again 
consider the purely electric Hamiltonian H + = ^-(E 2 + E 2 ), which commutes with 
the unitary operators 


T = exp 



+ —v 



T = exp 



-K% - — f 


9 

2n 



(2n i ~ i~\ 

= exp \ T dv + P - f )' 

= exp {^-df - j p - J , > 'j , (5.25) 
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that induce translations by ^ of the wave function (up to an abstract gauge trans¬ 
formation), i.e. 


T^{ip, V?) = exp &(<p + y, cp), 

TV(<p, <p) = exp V(<p, <p + y ). (5.26) 


It is straightforward to rewrite the boundary condition eq. (15. Ill) as 




(5.27) 


This immediately implies that T and T respect the domain structure of the Hamil¬ 
tonian, i.e., after acting with these operators, the wave function still obeys the 
boundary condition eq. (15.111) that defines the domain of H + . The operators T and 
T obey 


TT = exp 



(5.28) 


These two operators generate a discrete group Q, known as the magnetic translation 


group 


50, 511, which consists of the elements 


g{n, n, m) = exp 


2'Kim\ 

nr) 


rj-rnrj-ir 


7i, n, m G {0,1,..., k — 1}. 


(5.29) 


The group multiplication rule is given by 

g(n,n,m)g{n',n' ,m') — g(n + n' ,n + n' ,m + m' — nn'). (5.30) 

Here all summations are understood modulo k. The unit element is given by 


1 = 3 ( 0 , 0 , 0 ), 


(5.31) 


and the elements 

z m = 0(0,0, m) = exp 

form a cyclic Abelian subgroup Z(fc) C Q. 
g{n,n,m) is given by 


(T) • « 5 -> 

The inverse of the group element 


g(n, 7i, m) 1 = g(—n, —n, —m — nn). 

This follows because 


(5.33) 


g(n, n, m)g(—n, —7 1 , —m — nn) = g( 0, 0, — nn + nn) = g( 0, 0, 0) = 1. (5.34) 

Let us now consider the conjugacy class of a group element g(n,n,m) that consists 
of the elements 


g[n!, n'■> m ')g( n , n, 7n )d( n \ n', fn') 1 = 

g{nl + n, n' + n, m' + 7 n — n'n)g{—n', —nl, —ml — n'n') = 

g(n, n,m — n'(n + n') + ( n! + n)n') = g{n, n,m + nn' — n'n). (5.35) 
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The elements <7(0,0, m) = z m G Z (k) are conjugate only to themselves and thus 
form k single-element conjugacy classes. Multiplication by a phase z m amounts to 
a Z {k) gauge transformation. The conjugacy classes hence correspond to gauge 
equivalence classes. The elements <7(0, 0,m) = z m commute with all other elements 
and thus form the center Z (k) of the group Q. Since the individual elements of the 
center form separate conjugacy classes, the center forms a normal subgroup which 
can be factored out. This corresponds to identifying gauge equivalence classes. It 
should be pointed out that Q is not the direct product Z (k) x Z (k) x Z (k). Since the 
quotient space Q /Z(fc) is not a subgroup of Q , Q is not even a semi-direct product 
of Z {k) x Z (k) and Z (k), but just a particular central extension of Z {k) x Z {k) by 
the center subgroup Z {k). 


5.5 Spectrum of the Hamiltonian on a Single Cross 


Before we consider the full Hamiltonian of the entire lattice theory, let us again 
consider the purely electric Hamiltonian H + of a single cross. It commutes with 
both T and T, which, however, don’t commute with each other. Let us construct 
simultaneous eigenstates of H + and T. Since T k — 1, the eigenvalues of T are given 
by exp(27u///c) with l e {0,1,...,A; — 1}, while the eigenvalues of H + are given 


by E n = M(n + |) (where n e {0,1, 2,...} denotes the Landau level 52, [5 


m 


the “mechanical” analog, and M is the “photon” mass). Hence, we can construct 
simultaneous eigenstates \nl) such that 


H + \nl) 


M 



nl ), T\nl) = exp 



nl). 


(5.36) 


As a consequence of eq. (l5.28p one obtains 


TT\nl) = exp 



TT\nl) = exp 


i(l — 1) ^ 


T\nl), 


(5.37) 


which implies 

T\nl) = | n(l — 1)). (5.38) 

Since [T, H] = 0, the k states | nl) with l e 0,1,..., k — 1 are degenerate and form an 
irreducible representation of the magnetic translation group. If we restrict ourselves 
to the lowest Landau level, n = 0, by sending e — » 00 such that M —» 00 , the local 
Hilbert space of a cross is Axdimensional. As we will see later, this is the limit in 
which the doubled Chern-Simons-Maxwell theory reduces to the Z (k) variant of the 
toric code. 


In order to prepare for the derivation of the toric code in a later subsection, it 
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is useful to explicitly construct the eigenfunctions | nl), which are given by 


V) = {<P<p\nl) = -7= Y ^ ( V ~ 27r 


y/2n 

x exp I i{ktp + 9) 


mEZ 

l 9 
m + — + 


l 9 
m + — + 


k 2 tt k 


k 2irk 


kipip 

i —:— 

47T 


(5.39) 


Here 


w?) = V^y^ exp r^ lff ” 



2vr 


<P 


(5.40) 


(with H n denoting the Hermite polynomials) are the eigenfunctions of a 1-dimensio- 
nal harmonic oscillator with characteristic “frequency” given by the “photon” mass 
M = The corresponding energy eigenvalue is thus given by E n = M(n + |). 
It is straightforward to convince oneself that these wave functions indeed satisfy 
the correct boundary conditions eq. 115. lip and are eigenstates of H + . The boundary 
conditions on the torus parametrized by the self-adjoint extension parameters 6 and 
9 are given by 


£) = £) = (fi). (5.41) 

The boundary conditions of eq. 115.411) take the form of’t Hooft’s U(k) twisted bound¬ 
ary conditions [54, [55|. The transition functions 11,(1 6 U(k) play the role of U(k) 
gauge transformations. Their matrix elements are given by 


kli V = exp 


.27 rl + 9 

~k~ 


hi', ^u 1 — di+i,i', 


(5.42) 


where di+iji is understood modulo k. The transition functions satisfy the cocycle 
consistency condition 


fin = exp 



nn, 


(5.43) 


which is characterized by the U(k) gauge invariant twist factor exp(2m/k) G Z (k). 
The twist factor is an element of the center of U(k). The operators fl and H obey 
the same commutation relation as T and T (cf. eq. ll5.28p ) and thus also generate 
a representation of the magnetic translation group. As 9 is increased by 27r, the 
Z (k) electric flux quantum number l increases by 1. An adiabatic change of the 
self-adjoint extension parameters thus allows us to continuously interpolate between 
the different discrete Z (k) electric flux sectors. 
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5.6 U(k) Berry Gauge Fields 


It is natural to ask how the wave functions on a single cross respond to an adiabatic 
change of 9 and 9. Since k states \nl) (with l e {0,1,..., k — 1}) are degenerate, 
this gives rise to an abstract U(k) non-Abelian Berry gauge held on the torus f/(l) 2 
parametrized by 9 and 9 , 

9w(0, 0) = (nl\W ] d e W\nl') = 0, W = exp ^(#y + 0<p)j , 

gw(9,9) = (nl\W^dgW\7il') = i^-6 w . (5.44) 

In order to account for the fact that different values of 9 and 6 correspond to different 
domains of the Hamiltonian, it is important to include the unitary transformation 
W in the definition of the Berry connection. The gw and 'gw are elements of two 
k x k anti-Hermitean matrix-valued fields g and g that play the role of abstract U(k) 
Berry vector potentials. As we see, g turns out to be Abelian, and g even vanishes, 
at least in the gauge that we have picked. The corresponding Berry held strength, 
i.e. the abstract “magnetic” held, is given by 

MM) = *j(M) - &ss(0,0) + 19(9,9),9(9,9)], MM) = (5.45) 

It is again Abelian and constant over the 2-dimensional torus parametrized by 9 and 
9. Since the held strength is proportional to the unit-matrix, it is even invariant un¬ 
der general non-Abelian U(k) gauge transformations (resulting from a basis change 
in the subspace of degenerate states | nl)). This means that it actually reduces to an 
Abelian 17(1) Berry held strength. Still, since the held strength is non-zero, it can 
give rise to non-trivial Berry phases. However, since it is Abelian, its use for quan¬ 
tum information processing will be limited. The total “magnetic” hux that threads 
the torus is given by This is consistent with the twist factor exp(2iri/k) e Z (k) 
of eq. (15.431) that characterizes’t Hooft’s twisted boundary condition [54, 155!] • 


Besides the Berry held strength, there are gauge invariant Berry Polyakov loops 
that result from the holonomies of the torus parametrized by the periodic self-adjoint 
extension parameters (9,9). By integrating the Berry connection along a straight 
line wrapping around the torus, we obtain 


-2 IT 


d9 g u (9,9) = 0, 


p2tt _ ^ n 

d9 gu(9,9) = i 


k 


(5.46) 


The Polyakov loop matrices <$>,<$> E U(k), which describe parallel transport around 
the torus, also receive contributions from the transition functions and H such that 


®w(0) = V,iv = exp 


.27rZ + 9 

~k~ 


ill 1 , 


®w(9) = exp 





exp 



Si+1,1'- 

(5.47) 
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As we will see later, this will allow us to derive the mutual statistics angle ^ also 
in the compact Chern-Simons theory. 


5.7 Z (k) Gauge Symmetry and Gauss Law 


The operators T and T are associated with the links of a given cross. Until now 
we have considered a cross centered at x £ X. Now we also add the corresponding 
objects for x G X such that 


T X}i = exp 
T Xji = exp 


27t i _ i ~ . 

+ 2 €i ^ x - J ~ J. €i P x -J ) > 


2tt i i 

+ ~r e v^x,j 


(5.48) 


A Z (k) Gauss law combines T or T operators from links that touch a common lattice 
point to form unitary operators that represent local Z (k) gauge transformations 


V x = 


IF 


d 


x+2 M 


X 


G A, 


V x = 


U f 


x+-=i,i 


X 


G A. 


(5.49) 


By construction, V x and V x represent the manifest Z (k) gauge symmetries of the full 
Hamiltonian (including both the electric cross and the magnetic plaquette contribu¬ 
tions), which respect the domain structure defined by the values of the self-adjoint 
extension parameters 6 and 6. Although T and T on the same cross do not commute, 
original and dual gauge transformations on neighboring sites commute (cf. Fig J3]) 


LrL J , + aj + |2 — 

T 1 T Tt T 't T T T 1 ! T 1 ! _ 

1 x+%l,l 1 x+%2,2 1 x -!ii, 1 1 x-p,2 i *+ol+f2,l 1 x+%l+a2,2 1 x+ !t%l 1 x+%1,2 ~~ 

T ■ - T - - T - T - — V - ~V 

J x+al+f2,l J a; +fl+a2,2 J a;+ |2,l J 3; +fi,2 J ^+fl,l J ^+f2,2 J a ._|I il J 3; _|2,2 _ G+fl+f2^' 

(5.50) 


Here we have used 


T x+%l,l T , 


T 


4 

x+f i,2 

it 


= exp —— 


;+«2.2U+I2.1 = exp 


—j ft T . 

2vri \ ~ t 


jfc 


T t t 

1 x+ 0. 1 1 x+p,2- 


(5.51) 


For given 6 and 9, physical states without external charges from the corresponding 
domain of H must obey the Z(k) Gauss law 

V X \V) = |tt), x G A, V X \V) = |T), xG A. (5.52) 
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States that have external charges Q X ,Q X £ 0,1,..., k — 1 on the original and dual 
lattice obey 



(5.53) 


As we have seen before, 6 and 6 represent external non-dynamical U( 1) lattice gauge 
fields. When we perform U( 1) gauge transformations, we may change the values 
of 9 and 9 , which leads us into a unitarily equivalent domain of H , thus leaving 
the physics invariant. The physics does, however, depend on the gauge invariant 
plaquette field strengths exp (irj x ) and exp {irj x ) of eq. 05.24|h 

5.8 Pure Chern-Simons Limit 

Let us now take the limit of large e 2 , which implies that the “photon” of mass 
M = is removed from the spectrum and the system reduces to a pure Chern- 
Simons theory. In particular, the low-energy physics is restricted to the lowest 
Landau level n — 0. Then there are only k states left on each individual cross. As 
before, we choose an electric flux basis of the local Hilbert space, characterized by the 
fluxes l £ {0,1,..., k — 1} on the links of the original lattice. The electric field energy 
is then given by for each cross, and the dynamics is entirely controlled by the 
magnetic field energy as well as by the Gauss law. Since they are suppressed by 
the magnetic contributions to the energy can then be treated in first order degenerate 
perturbation theory, by evaluating their matrix elements between the degenerate 
electric flux eigenstates |0Z). Putting n — 0 in eq. (15.391) . it is straightforward to 
obtain 


(± <V ±1 = Cexp (0/|T Tl |0O, 


(0/| exp(±f(^)|0/') = Cexp(±— 



(5.54) 


Using T|0Z) = 10(Z — 1)), this implies that (up to a trivial additive constant) the 
effective low-energy Hamiltonian (in the n — 0 sector) is given by 



(5.55) 
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where 77 and rj are the plaquette field strength variables associated with the non- 
dynamical gauge fields 6 and 9 of eq. fl5.24p . The first contribution results from a 
cos (a 2 B x ) term on the original plaquettes and the second contribution originates 
from a dual plaquette term cos (a 2 B x ). When the Hamiltonian acts on eigenstates 
| Q, Q) with charges Q x , Q x on the original and on the dual lattice (cf. eq. (l5.53p ) one 
obtains 


H e s\Q,Q) = 


C 4 


e 2 a 2 


xGA 


cos 



I Q,Q). 


(5.56) 

Hence, the effective Hamiltonian punishes violations of the Z(/c) Gauss law on the 
original and on the dual lattice. In particular, for rj x — rj x — 0, the ground state 
is reached when I4|T) = |T) for x G A and I4| v b) = |'F) for x G A, i.e. when 
Qx — Qx — 0 . For other values of rj x and rj x other constellations of charges are 
energetically favored. 


Let us construct the states | Q,Q) more explicitly. For this purpose, we construct 
dual charge projection operators at each site x G A 


1 


k -1 


Qx ^ ^ 

?n =0 


exp 


2nim 


k 


which obey 




2ni 


P Qx P Q' x ~ ^Qx,Q' x P Qx' VxP Qx — ex P ( Q* ) P Qx 


(5.57) 


(5.58) 


We now consider a state |[Z]) corresponding to a specific configuration [/] of electric 
fluxes that obeys the Gauss law in the presence of charges Q x G Z (k) on the original 
lattice, igA, i.e. 

~ l x-fi,i) modk = Qv ( 5 - 59 ) 

i 

By construction, we then have 

V x \[l])=exp^^Q^\[l]), x G A. (5.60) 


Next we apply the charge projection operators at all dual sites to obtain 

IQ.Q)=vn/> s J[i]), ( 5 . 61 ) 

xgA 

where AT is a normalization factor. By construction, this state indeed obeys 

V X \Q, Q) = exp (-jrQx'j IQ, Q), X e A. (5.62) 
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As we have seen in FigJU in the non-compact theory the transport of a charge 
along a curve C = dS on the original lattice, that encircles a total charge Qs = 
^2 xgsc a Qx 011 the dual lattice, is associated with a non-trivial topological phase 
exp(27 riQs/k), thus showing that original and dual charges have mutually anyonic 
statistics with the statistics angle ^ 7 . We will now derive the same result in the 
compact theory. For this purpose we again consider a Wilson loop 

Wc — n u x,i = n ex P(5.63) 
(. x,i)eC ( x,i)ec 

Using eq. (15.541) we then obtain 

Wc\Q,Q) = C L < V* II ew^)\Q,Q)=C Lc exp(^^ + i^\\Q,Q). 

x^ScA {x,i)GC \ / 

(5.64) 

Here Lq is the length of the closed loop C and Qs = J2 X £S<zaQ x the total dual 
charge encircled by C. As before, the phase of the Wilson loop contains the Q 5 -term 
which implies that original and dual charges have mutually anyonic statistics with 
the statistics angle Due to the non-trivial background gauge held 9, there is an 
additional phase determined by the total background flux i]s = X^eScA % encircled 
byC. 

The Wilson loop describes the instantaneous transport of a charge around the 
loop C. Alternatively, we now want to consider a much slower adiabatic process in 
which a charge-anti-charge pair is created, transported around the loop C, and finally 
annihilated. For this purpose, we compute the Berry phase that we accumulate when 
we gradually change the background held 6 from 0 to 2tt on all dual links intersecting 
the loop C. According to eq. (15.471) this provides us with a factor exp (i9 x>i /k) for all 
links that belong to C. By Stokes theorem, this again yields the phase exp (irjs/k). 
In addition, the factor 5i + \^ in eq. (15.471) shifts all electric hux variables along C 
by 1 , which again gives rise to the phase exp(27ri(5s/fc) when one returns to the 
initial state | Q, Q) at the end of the adiabatic charge transport process. This shows 
explicitly that the anyon statistics angle ^ is insensitive to the details of how the 
charge transport process is realized, be it adiabatic or instantaneous. 


5.9 Relation to the Toric Code 


The toric code is a Z(2) lattice gauge theory with a degenerate ground state that can 
be used as a storage device for quantum information that is topologically protected 
against decoherence. The fundamental degrees of freedom on each link of a square 
lattice are quantum spins |. In fact, the toric code is a simple member of^ large 


class of unconventional lattice gauge theories known as quantum link models 56il60 
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The Hamiltonian of the Z(2) toric code is given by 


H 


xeA 


xeA 


V x = exp (nrS* exp 



(5.65) 


There is a spin |, S x . tl associated with each link connecting the neighboring lattice 
sites x — § i and x+%i. The quantum link operator is given by U Xji = S xi , such that 
the plaquette term 


U x = Si 


x-p,l 


nl nl nl 

z+fi,2 z+§2,l ~ 


x-fl,2’ 


x e A, 


(5.66) 


is invariant under discrete Z(2) gauge transformations kl x = ±1. The G'-term with 
G > 0 punishes violations of the Z(2) Gauss law V x |T) = IT). Note that for Z(2) 
one has U x = U\ and V x = Vj. 


One can describe the dynamics of the toric code in a basis of electric flux states 
on each link, which are characterized by the 3-component j = ±| of the spin. The 
two spin values then correspond to two directions of the electric flux on each link. 
The quantum link operator U Xti = S x i reverses the direction of the electric flux, and 
thus the plaquette operator U x reverses the direction of all electric fluxes encircling 
the plaquette. The Z(2) toric code naturally generalizes to a Z (k) variant with k 
distinct values of the electric flux on each link. The corresponding quantum link 
operator then increases the flux by one unit modulo k. The Z (k) Gauss law ensures 
that the fluxes entering a common lattice point add up to zero (again modulo k). 
Interestingly, the Z (k) variant of the toric code is exactly what emerged from the 
doubled compact Chern-Simons-Maxwell theory in the pure Chern-Simons limit of 
large “photon” mass, when we identify J = G = 


5.10 Rotating 9 and 6 into the Hamiltonian 

Until now, the external background gauge fields 9 and 9 entered the theory via 
the boundary conditions eq. (15.111) that define the domain of the Hamiltonian. For 
completeness we would now like to show that the background fields can be rotated 
into the Hamiltonian by the unitary transformation 

w = n exp ( —^{9 x ,l<Px ,2 + 0 X ,2<Px,l) ) II exp ( + Ss.l^x^) J , 

H' = WHW\ |T') = VF|T). (5.67) 
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The original Hamiltonian 

„2 


H = 


e 

~2 


E 

xex 

+ tE 






e 2 a 2 


^ cos {a 2 B x ) 

xeA 


e 2 a 2 


^cos (a 2 B x ), 


(5.68) 


iGA 


with the boundary condition (for a cross with x E X) 

k .. ~\ 


T(<^ + 27t, <p) = exp ( —+ id ) T((p, <p), T((p, <p + 27r) = exp ( T(<^, <p), 


2 y 

then turns into the new Hamiltonian 
H' = 


k 


(5.69) 


e 

~2 


+ 


e 

~2 


E 

xex 

E 

xex 


■ q k __ 0 X ,2 

+ 27 


k ~ #2-1 

- s9 ^,> + 4j^.i - 


+ (~ i9 « + + f^) 

+ (~ iS « _ _ If) 


^ E cos(a 2 B„) - -h E ™s(a 2 BJ, 

a^GA ;ceA 


(5.70) 


The new Hamiltonian depends on 0 and 0, while the new wave function now obeys 
the simplified boundary condition (again for a cross with x E X) 

’I'' [ 'p + 2ir, ip) = exp (-4p) <p'(p, p). 'J'V. p + 2?r) = exp (4p) <p'(p, pj. 

(5.71) 

which is 0- and d-independent. The wave function on a single cross x E X then 
takes the form 


^'nf (<A <p) = {<p<p\nl) = —= Y V'n ( <P - 2vr 


mEZ 



l 

6 

1 

to 

m + — + 

k 

2nk 


x exp z(£> 


km + / 


kip 
47r 


+ i 


/ 0 
m + — + 




fc 27 tA; 27t 


0 . 


(5.72) 

In this case, in contrast to eq. fl5.44p . and using |nZ)' = W|nZ), the definition of the 
Berry connection no longer requires the insertion of W and is simply given by 

9w(d, 0) = '(nl\d e \nl'y, g w (9, 9) = '{nl\d~\nl')'. (5.73) 

The resulting values of the Berry connection are the same as in eq. (15.441) . 
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6 Conclusions 


We have obtained the Z(k) variant of the toric code in the infinite “photon” mass 
limit of a doubled U( 1) Chern-Simons-Maxwell lattice gauge theory. In contrast 
to ordinary lattice gauge theories, in this case the held algebra is not link-based. 
Instead it is based on a cross formed by a link and its corresponding dual link. In 
ordinary lattice gauge theory, each individual link has a “mechanical” analog: a 
“particle” moving in the group space, i.e. a “particle” moving on a circle for a com¬ 
pact Abelian 1/(1) lattice gauge theory. In the cross-based doubled compact Chern- 
Simons-Maxwell theory, on the other hand, the “mechanical” analog is a charged 
“particle” moving on a torus I/(l) 2 which is threaded by an abstract “magnetic” 
flux, determined by the prefactor k of the Chern-Simons term. The Dirac quanti¬ 
zation condition for the abstract “magnetic” flux then leads to the quantization of 
the level k E Z. 

Remarkably, the cross-based electric contribution to the Hamiltonian has two 
self-adjoint extension parameters 9 and 9, which themselves form two non-dynamical 
background 1/(1) lattice gauge fields. In a fixed background, the manifest gauge 
symmetry of the compact Chern-Simons-Maxwell theory is then reduced from I/( 1) 
to Z {k). In particular, only a Z (k) (and not the full U{ 1)) Gauss law must be im¬ 
posed on physical states in a given domain of the Hamiltonian (which is determined 
by 9 and 9). The explicit symmetry breaking from U( 1) to Z [k) has a quantum 
mechanical origin, because only the quantum but not the classical theory is sensi¬ 
tive to the external self-adjoint extension parameters, which manifest themselves as 
Aharonov-Bohm phases. One might thus describe the quantum mechanical gauge 
symmetry breaking as an “anomaly”. However, it is more like the explicit breaking 
of CP invariance due to a non-zero 0-vacuum angle in 4-dimensional non-Abelian 
gauge theories. It is remarkable that a gauge symmetry can be broken by quantum 
effects in a similar manner. It should, however, be pointed out that the redundancy 
associated with the 1/(1) gauge symmetry of the classical theory still persists at 
the quantum level. It manifests itself in the 1/(1) gauge redundancy of the exter¬ 
nal background lattice gauge fields formed by the self-adjoint extension parameters 
9 and 9. In this paper, we have not attempted to take the continuum limit by 
sending e 2 a —> 0. It is conceivable that the full U( 1) gauge symmetry would then 
emerge dynamically from the manifest Z (k) symmetry, and that the lattice theory 
would reduce to the non-compact doubled Chern-Simons-Maxwell continuum theory 
discussed in Section 3. 

While the electric held energy does not depend on the values of the self-adjoint 
extension parameters, the magnetic held energy does. The Z(/c) toric code emerges 
from the doubled Chern-Simons-Maxwell theory in the limit of infinite “photon” 
mass, M — y~ —> oo. In this limit the magnetic held energy vanishes, and thus 
the dependence of the energy on the self-adjoint extension parameters disappears. 
However, the wave function still depends on 9 and 9. As a result, the toric code 
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has a large variety of hidden self-adjoint extension parameters, which form two U( 1) 
gauge fields, one associated with the original and one associated with the dual lattice. 
Under adiabatic changes of the parameters 9 and 9, a U(k) Berry gauge field with 
non-trivial Berry curvature arises. Since the Berry field strength turned out to be 
gauge equivalent to an Abelian field strength, manipulating the corresponding Berry 
phases will certainly not allow universal topological quantum computation. Still, it 
will be interesting to investigate further whether manipulating these parameters can 
be utilized for other forms of quantum information processing. 


In this work, we have encountered a large variety of different gauge fields. First 
of all, we started out with non-compact Abelian gauge fields A and A, both in the 
continuum and on the lattice. Upon compactification of the gauge group, this gave 
rise to the compact 17(1) lattice gauge fields <p and <p. The group space of an original 
and dual cross-based link-pair f/(l) 2 was endowed with an abstract vector potential 
(a, a) describing the quantized “magnetic” flux represented by the level fceZ. The 
Polyakov loops winding around the group space torus U( l) 2 then gave rise to two 
self-adjoint extension parameters 9 and 9 , which play the role of Aharonov-Bohm 
phases for the “mechanical” analog “particle” moving in the group space. While 
they are invisible at the classical level, at the quantum level the parameters 9 and 
9 break the manifest gauge symmetry from 1/(1) down to Z(k). Still, the full 17(1) 
redundancy of the gauge description remains, because the self-adjoint extension 
parameters themselves turned out to be non-dynamical background lattice gauge 
fields associated with the original and dual lattice, respectively. An extended version 
of the toric code then emerges at low energies, which is a Z (k) lattice gauge theory. 
Finally, by considering adiabatic changes of the self-adjoint extension parameters 
9 and 9, U(k ) Berry gauge fields arise, which are defined over the base space of 
the lattice gauge fields 9 and 9. All these different gauge structures are intimately 
related to one another, and deserve further study, in particular, in the context of 
quantum information processing. 


We have often mentioned the “mechanical” analog of a “particle” moving in the 
cross-based group space torus U( l) 2 , threaded by k units of quantized “magnetic” 
flux. It would be most interesting to build such a system in the laboratory, in 
order to manipulate quantum information. While building a torus seems like a most 
difficult task (at least to us), we have great confidence in the ingenuity of AMO 
experimentalists. Perhaps one can manipulate trapped ions to mimic the physics 
of the cross-based “mechanical” analog. Using digital quantum simulations with up 
to 100 gate quantum operations (oij |. this has already been achieved for a single 


plaquette of the ordinary (i.e. not extended) Z(2) toric code 62]. If one could 


experimentally incorporate the background gauge fields 9 and 9 , putting together 
several crosses would provide us with a remarkable quantum system that embodies a 
large variety of Abelian dynamical or background gauge fields as well as an additional 
U(k) Berry connection. 
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It would be interesting to extend our investigations to non-Abelian Chern-Si- 
mons-Maxwell theories on the lattice. In particular, one may ask whether they also 
reduce to variants of the toric code with a discrete non-Abelian gauge group. We 
leave this problem for future studies. 
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